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Linear magnetoresistance due to multiple- 
electron scattering by low-mobility islands 
in an inhomogeneous conductor 
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Linear transverse magnetoresistance is commonly observed in many material systems 
including semimetals, narrow band-gap semiconductors, multi-layer graphene and topological 
insulators. It can originate in an inhomogeneous conductor from distortions in the current 
paths induced by macroscopic spatial fluctuations in the carrier mobility and it has been 
explained using a phenomenological semiclassical random resistor network model. However, 
the link between the linear magnetoresistance and the microscopic nature of the electron 
dynamics remains unknown. Here we demonstrate how the linear magnetoresistance arises 
from the stochastic behaviour of the electronic cycloidal trajectories around low-mobility 
islands in high-mobility inhomogeneous conductors and that this process is only weakly 
affected by the applied electric field strength. Also, we establish a quantitative link between 
the island morphology and the strength of linear magnetoresistance of relevance for future 
applications. 
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Linear transverse magnetoresistance (MR) is a phenomenon 
commonly observed in many material systems, including 
semimetals 1 , narrowband-gap semiconductors 2,3 , multi-layer 
graphene 4 and topological insulators 5 . Linear MR can arise from 
either quantum 6 or classical effects 7 ' 8 . In particular, it can originate 
in an inhomogeneous conductor from distortions in the current 
paths induced by macroscopic spatial fluctuations in the carrier 
mobility; this type of non-saturating linear MR has been observed 
in several materials and explained in terms of a phenomeno- 
logical semiclassical random resistor network model 7 ' 8 . However, 
the link between the linear MR and the microscopic nature of 
the electron dynamics remains still unexplored. Existing theories 
of classical linear MR apply to systems with ohmic conductivity, 
whereas the effects of hot electron dynamics on the linear MR are 
largely unknown. On the other hand, a significant enhancement 
of MR has been observed in the non-omhic space-charge limited 
transport regime 9 and under avalanche breakdown 10 , thus rais- 
ing problems of fundamental interest and technological relevance, 
that is, what is the specific carrier dynamics giving rise to the linear 
MR, is it preserved in the presence of a large applied electric field 
and can it be controlled and exploited in sensitive Hall sensors or 
magnetoresistors 1 1 ? 

Here, we examine the physical origin of the classical linear MR 
in a high-mobility inhomogeneous conductor. Our measurements 
and Monte Carlo simulations focus on the n-type compound semi- 
conductor InAs and consider an extended range of applied electric 
(E) and magnetic (B) fields, and different sample geometries. By 
modelling our data with Monte Carlo simulations of the electron 
dynamics, we demonstrate that the linear MR arises from multiple 
scattering of the current- carrying electrons by low-mobility islands 
within the conducting layer. This process is only weakly affected by 
the details of the electronic band structure and electric field strength. 
Also, it persists at magnetic fields as high as -50 T corresponding to 
a magnetic length as small as 4nm, significantly smaller than the 
electron mean free path (-100 nm) and island size. The link between 
the linear MR and the stochastic behaviour of the electronic cycloi- 
dal trajectories around low-mobility islands is a finding of general 
interest as it is relevant to all high-mobility inhomogeneous con- 
ductors. The identification of a quantitative link between the island 
morphology and the strength of linear MR is also potentially useful 
for tailoring the linear MR in future applications. 

Results 

Material systems and devices. Our InAs and In(AsN) epilayers were 
grown by molecular beam epitaxy (MBE) on a semi-insulating (100)- 
oriented GaAs substrate. Owing to the mismatch between the substrate 
and epilayer lattices, inhomogeneities can arise due to the presence of 
dislocations and defects; also, the N-atoms act as scattering centres for 
the conduction electrons 12 ' 13 . Thus with increasing N-content from 0 
to 1 %, the electron Hall mobility decreases from /i H = 3.3m 2 V _1 s _1 
(3.1m 2 V- 1 s- 1 ) to 0.6m 2 V- 1 s- 1 (0.5m 2 V- 1 s- 1 ) at T=300K 
(4.2 K). Correspondingly, the electron density increases from 
rc e = 4xl0 16 cm~ 3 in InAs to n e = 1.4xl0 17 cm~ 3 in In(AsN), with 
both values being very weakly affected by temperature. 

The presence of inhomogeneous disorder combined with the 
high electron mobility makes these InAs and In(AsN) layers par- 
ticularly well suited to the study of classical linear MR: even in the 
low-mobility layers the high values of jlB > 1, which are necessary 
for observing a large MR, can be achieved with available conven- 
tional (superconducting) solenoids. To investigate the influence of 
sample geometry and electric field on the MR, two -terminal planar 
devices with InAs channels of length L = 2, 5 and 10 |im and of width 
W= 1, 2 and 5 |im were also studied. 

High magnetic field studies. In our experiments, the magnetic 
field was applied parallel to the growth axis z, that is, B= [0, 0, B z ]. 
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Figure 1 1 Linear transverse MR in Hall bars and short-channel devices. 

(a) Transverse MR, Ap xx /p xx , as a function of magnetic field, B ZI for an 
InAs Hall-bar at 7 = 4.2, 100 and 300 K. (b) Ap xx /p xx as a function of B z 
for InAs (black) and In(AsN) (green) Hall-bars (T= 300 K). Dashed lines 
and arrows illustrate the transition of the dependence of Ap xx /p xx on B z 
from quadratic to linear at the characteristic magnetic field B z l . A linear 
fit to the MR data by the relation Ap xx /p xx = iiqB z (s/2I) gives an average 
island coverage factor f~0.M and 0.32 in InAs and In(AsN), respectively. 
Here, s is the average size of the low-mobility regions, / is the inter-island 
separation and {Iq is the zero field mobility, (c) Ap xx /p xx versus B z for short 
InAs channels of length, i, and width, W, at T= 100 K (line) and 300 K 
(symbols). The insets in panels (a) and (c) sketch a Hall-bar and a two- 
terminal planar device, respectively. 



The transverse MR ratio defined as Ap xx /p xx = [p xx (B z ) -p X x(®)]/pxx(Q) 
was studied with magnetic fields, B z , up to 14T, generated by 
a superconducting magnet, or in a pulsed solenoid at B z up 
to 47 T with a total pulse length of 100 ms. Figure la shows the 
^-dependence of Ap xx /p xx for an Hall-bar based on n-type InAs. 
A positive linear MR is observed at B z > 1 T and is weakly affected 
by temperature in the range T= 4.2-300 K. A weak oscillatory com- 
ponent of p xx is observed at low temperature (T<20K), which we 
attribute to the Shubnikov-de Haas effect in the two-dimensional 
electron gas that accumulates on the surface of InAs due to the bend- 
ing of the conduction band by charged native surface defects 12 . The 
resistance of the surface two-dimensional electron gas at B z = 0T is 
at least one order of magnitude higher than p xx and, at all values of 
B z and T, the amplitude of the magneto-oscillations is much smaller 
than the measured monotonically increasing MR. 

Linear MR can occur when the electron distribution enters the 
so-called quantum limit, in which only the lowest Landau level is 
occupied, that is, when the cyclotron energy, hco c , exceeds the Fermi 
energy, £ F , and when hco c , £ F >A: B T (ref. 6). As shown in Fig. lb, 
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Figure 2 | Monte Carlo simulations of electron trajectories in magnetic field. Calculated electron trajectories (red lines) in InAs for an applied electric 
field E x = 0.5kVcm _1 and under various applied magnetic fields B z : (a) B z = 0, (b) 1, (c) 10 and (d) 20T. The black areas correspond to regions of low 
mobility. As shown in the magnified trajectories in the insets of panels (b) and (d), in the vicinity of a low-mobility island, the electron undergoes multiple 
scattering events, which act to deflect the cycloidal motion. The bar in panel (a) corresponds to 1 um. 



the transition of the MR to a linear B- dependence occurs at B z ~l T 
below the value satisfying the quantum limit condition hco c = £ F , 
that is, B z =h(37i 2 n e ) /3 /2e. This limit corresponds to fields >3T 
for an electron density n e = 4xl0 16 cm -3 , which we derive from 
the gradient of the Hall resistivity p xy) in the limit B z -^0. Also, we 
note that the condition for single occupancy of the lowest Landau 
level requires that k^T<^hco c . This condition is satisfied for 
temperatures < 50 K at B z = 1 T and for the electron cyclotron mass 
of InAs (m e = (0.025±0.01)m 0 , where m 0 is the electron mass in 
vacuum) 14 . However, the linear MR is observed well above T= 50 K 
(see Fig. 1), indicating that it does not have a quantum origin, but 
is a classical effect. 

Classical linear MR has been explained using a phenomeno- 
logical model in which the distorted current paths arising from 
disorder-induced ihnomogeneities and macroscopic variations, 
AjU, in the electron mobility are simulated by a network of ran- 
dom resistors. Each resistor element has four terminals to take into 
account not only the local resistive voltage drop but also the Hall 
voltage associated with the local mobility, /I, of a particular region 
of the sample 7 ' 8 . The crossover of the MR from quadratic to linear 
Independence is predicted to occur at a characteristic magnetic 
field B z = JT l for AfiljK 1, or at B z f = Aji~ l for Afilji> 1, where jl 
corresponds to an average of the spatially varying mobility. A com- 
parison of the B z - dependence of Ap xx /p xx in our InAs and In(AsN) 
Hall-bars at T=300K (see Fig. lb) shows that both the value of 
II B z and the slope of the MR follow the systematic decrease of jl 
from 3.3m 2 V" 1 s" 1 in InAs to 0.6m 2 V" 1 s" 1 in In(AsN). The 
characteristic magnetic field B z provides an estimate of the aver- 
age mobilities jl = \IB Z ~ lm 2 V~ 1 s~ 1 and 0.5m 2 V _1 s _1 in the 
InAs and In(AsN) samples, respectively, which are comparable with 
those measured by the Hall effect. Also, the MR tends to follow the 
temperature dependence of jl: in the InAs sample, an increase of 
T>4.2K and up to T -100 K leads to an increase of /i H from 3.1 to 
3.8m 2 V _1 s _1 followed by a decrease at T>150K. This variation 
is caused by the interplay between impurity and phonon scatter- 
ing and is accompanied by a similar T- dependence of the MR (see 
Fig. la). Linear MR is also observed in the two -terminal small 
channel devices, see Fig. lc. 

Monte Carlo simulations. As our InAs and In(AsN) epilayers are 
grown on highly lattice -mismatched GaAs, threading dislocations 
tend to form at the epilayer/substrate interface. The macroscopic 
(>0.1jLim) inhomogeneities observed by transmission electron 
microscopy of our samples 15 are likely to induce variations in the 
electronic properties. We model the effect of such inhomogenities 



on the electron dynamics by Monte Carlo simulations 16 . As the mag- 
netic field is applied along the z-direction, the magnetic field affects 
the components of the electron motion in the x-y plane. In our sim- 
ulation of the two-dimensional electron trajectories, inhomogeneity 
is introduced by placing low-mobility islands randomly within the 
x-y plane (Fig. 2). The random spatial profile of the islands, A(x,y), 
is generated from the power spectrum of the autocorrelation func- 
tion (A(r)A(r')) = A 2 exp[- |r-r'| 2 /A 2 ], where A is the correlation 
length, r=(x, y) and r' = {x\ /). The low-mobility regions R\ ow 
are then defined according to the relation £ low ={ r\A(r) < A th }. 
This condition determines the fractional area,/, of the sample with 
low mobility in which electrons undergo scattering at a rate higher 
than in the remaining high-mobility part of the sample. In a three- 
dimensional model, the value of/ would result from an average of 
the /-values from various x-y planes of the three-dimensional con- 
ducting channel. In the low-mobility regions, we use a scattering 
rate r is j an( j = r- Y + r e , where r { = 1.5xl0 12 s _1 and r e = 2.5xl0 14 s _1 are 
the elastic and inelastic phonon scattering rates, respectively; else- 
where, the electron drift is mainly limited by inelastic scattering with 
r~r v The non-parabolicity of the electron energy dispersion of InAs 
is modelled as £(k)[l + ae(k)] = h 2 k 2 /2m e , where ra e = 0.026m 0 and 
a=2.2 eV~ l . This form of the energy dispersion occurs in numerous 
semiconductors, including the narrow gap InSb for which observa- 
tions of linear MR have been also reported 2 ' 3 . The electron veloc- 
ity, given by v(k) = h~ l V k e, and the wave vector, k, change with 
time under the action of the Lorentz force, F = -e[E + vxB] y where 
E= [E x , E r 0] and B = [0,0, B z ] 17 . The Hall field component, E y , is set 
so that they- component of the average electron velocity is zero, that 
is, v = [v x ,0,0], consistent with the current flow for our Hall bars. 
The main findings of our simulations are not changed if we consider 
a different geometry in which the Hall field is effectively shorted 
out', that is, Ey = 0. Thus, our analysis applies not only to Hall-bars 
but also to short channels and Corbino geometries, in which the 
Hall electric field is fully short-circuited' by the contact electrodes. 

Figure 2a-d illustrates the calculated electron trajectories at 
£ % = 0.5kVcm _1 and B z = 0, 1, 10 and 20 T. Here, we have assumed 
a correlation length A = 0.5|im,/=5% and low-mobility regions of 
average size s = 0.4|im and inter-island separation / = 8|im. It can 
be seen that in the vicinity of the low-mobility islands, the electron 
undergoes multiple scattering events, which act to deflect the cycloi- 
dal motion (see inset of Fig. 2b). This island-induced scattering is 
observed at low (Fig. 2b) and high (Fig. 2d) applied magnetic fields. 
It is absent in a homogenous system and is responsible for the linear 
magnetic field dependence of the inverse mobility component, 
jl xx ~ l , whose strength increases with/, as shown in Fig. 3a. 
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Figure 3 | Monte Carlo simulations of the electron mobility and its 
magnetic field dependence, (a) Dependence of the inverse of the mobility, 
jU xx -1 on magnetic field, B z , in InAs at various values of the island coverage 
factor f = 0.1 (black), f = 0.05 (red) and f = 0.02 (blue), (b) Calculated 
electron trajectories (red lines) in InAs at B Z = '\T, /f x = 0.5 kVcm" 1 and 
different electric field components E v (left: E„ = 0 and v v = 5x10 4 ms _1 ; 



right: 



and v v = 0). The black areas show regions of low 



mobility. A and B indicate two typical islands, (c) Dependence of fx xx 1 
on E x for f= 0.05. 



Discussion 

We interpret the results of the Monte Carlo simulations by consid- 
ering a simplified Drude model of the average electron velocity v x 
and average hopping distance travelled by the conduction electrons 
along the x- direction. In a uniform system, that is, with no islands, 
and for E y = 0 and vy£0, the average velocity along x is v x = AX/t, 
where AX =m e \E x \/(eB 2 ) is the typical hopping distance and T is 
the scattering time. In the presence of inhomogeneity, an elec- 
tron travels an additional distance AX'~s/2 around a low-mobility 
island. As the electron hits an island every time interval T is i an( j~//v H , 
where v H = E x /B z is the electron Hall velocity along y and / is the 
distance between the islands, the average electron velocity v x and 
mobility component jl xx increase as v x - AX/T+ AX7T is i an( j and 
jd xx = m e /(ze 2 B 2 ) + s/(2lB z ), respectively. Thus, if we compare 
this result with the case of an homogenous system for which 
ji xx = mJ(Te 2 B 2 ) and v x = m e \E x \/(TeB 2 ), we conclude that at 
high B z island scattering tends to facilitate electron drift along the 
x- direction, perpendicular to the cycloidal motion along y (see 
electron scattering around island A in Fig. 3b). 

For the geometry in which Ey = 0 and vy£0, the resistivity com- 
ponent p xx is given by p xx = E x /J x - Xlne]!^. Using the depend- 
ence of jd xx on B z at large fields, that is, \i xx -AX / /v ls \ ZIi ^ = s/{2lB z ) i 
we derive a linear MR, that is, Ap xx lp xx ~B z ti*, where jl* = }1q{s/2 I) 
is and effective mobility and {Iq is the zero field mobility. A similar 
conclusion can be drawn for the Hall geometry, that is, for £y£0 
and v y = 0. In this case the resistivity component p xx is given by 
Pxx ~ I 1 xx ^ ne (l^xx + Mxy)]> where jl xy = B z ~ l . Using this expression 
in the limits of high B z , we obtain also in this case that Ap xx /p xx 
~ B z ji*. Note, however, that in this geometry, the island scattering 
acts to decrease the average electron velocity, v x , by deflecting the 



electron motion along y (see electron scattering around island-B 
in Fig. 3b). 

In summary, the magnetic field dependence of the mobility 
component jl xx due to island scattering generates a characteristic 
linear dependence of the MR on magnetic field, which is the same 
in both geometries, that is, £y£0 and £^ = 0. This analysis identi- 
fies the physical origin of the random variation of the resistive and 
Hall effect components of each four-terminal resistive element of 
the network model 7 ' 8 : the randomness arises from the stochas- 
tic behaviour of the electronic cycloidal trajectories around the 
low-mobility islands. 

From a linear fit to the measured MR data for the large 
InAs Hall-bar of Fig. lb by the relation Ap xx /p xx = B z fi*, we obtain 
i u )f = 0.33m 2 V" 1 s" 1 . For i u 0 = i u H = 3.3m 2 V" 1 s" 1 , this corre- 
sponds to s/(2r)=jd*/jd 0 = 0A and an average island coverage fac- 
tor f=s/(s + l) 1-0.1 7. Similarly, for the In(AsN) Hall-bar of Fig. lb 
we find that s/(2 /) = l^/fl 0 = 0.23 and/= s/(s + Z)~0.32. Thus, in both 
cases the average island separation, /, is significantly larger than the 
average island size, s. Although we cannot deduce the exact value 
of / or s, we can infer a reasonable estimate from our MR experi- 
ments in small channel devices, see Fig. lc. Here, it can be seen 
that devices with small values of the channel length L (2-10 jam) 
and/or width W (1-5 |im) exhibit much smaller values of Ap xx /p xx 
and in the smallest devices the MR tends to saturate at high B z . The 
linear MR decreases systematically as the channel width W and 
length L decrease to values of the order of 1 jam. Thus, we infer that 
the dimensions of these small devices are comparable to the average 
island separation. 

The Monte Carlo simulations indicate that different values of 
the band non-parabolicity parameter a produce similar linear MR 
behaviour. Thus, the specific details of the electronic band struc- 
ture do not have a critical role in the observation of linear MR. Per- 
haps surprisingly, although ji xx ~ 1 tends to increase with increasing 
E x , its linear dependence on B z is found to be weakly affected by 
electric field (Fig. 3c) and only at low-magnetic fields ( < 1 T) and 
high E x does the value of ji xx ~ 1 tend to decrease with decreasing B z , 
thus leading to a weak negative MR. This is caused by the interplay 
between electron cyclotron motion and island-induced scattering, 
which tends to be increased by E x , that is, electrons hit a low-mobility 
islands more frequently at higher speeds. However, overall this 
effect tends to be weak and, at large B z , the deflection of the cycloi- 
dal motion by the low-mobility islands and corresponding linear 
MR are preserved. 

The weak dependence of the linear MR on applied electric field is 
confirmed by our MR experiments in the short channel InAs sam- 
ples. Figure 4a shows the dependence of the current density / on 
E x at low temperature (T=4.2K) and at various B z (up to 14 T) for 
an InAs channel of length L = 10jim and width W=2[im. Similar 
dependences of / on E x and B z were found in devices with differ- 
ent values of L and W. At £ Z = 0T, the J(E X ) curve is approximately 
linear at low bias. The ohmic behaviour is then followed by a slight 
sublinear bias dependence and, at larger electric fields, by a sharp 
rise of the current that is characteristic of systems switching from a 
state of low to high conductivity due to impact ionization and ava- 
lanche breakdown 15,18 . At all values of E X) the magnetic field acts 
to suppress the current and to shift the threshold electric field for 
impact ionization to significantly higher values. From Fig. 4b, it can 
be seen that the MR remains linear over an extended range of E x 
up to ~3kVcm _1 , just below the threshold electric field for impact 
ionization. In contrast, in the regime of avalanche breakdown the 
MR is strongly enhanced and tends to depart from a simple linear 
dependence due to the effect of the magnetic field in suppressing the 
avalanche multiplication of carriers 10 . 

To conclude, our data and analysis indicate that the presence of 
low-mobility islands in a high-mobility conductor and their effect 
on the electron trajectories underpin the observation of classical 



NATURE COMMUNICATIONS | 3:1097 | DOI: 10.1038/ncomms2106 | www.nature.com/naturecommunications 
© 2012 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2106 



ARTICLE 




e z (T) 



Figure 4 | Linear transverse MR at high electric fields, (a) Dependence 
of the current density J on electric field E x at various magnetic fields B z 
for an InAs channel of length /_ = 10|im and width \A/ = 2|im (7=4.2 K). 
Measurements of the J(E X ) characteristics were performed in both the 
direct current (lines) and pulsed mode (symbols). The inset sketches a 
two-terminal planar device, (b) Transverse MR, Ap xx /p xx , as a function of 
B z at various values of E x . The resistivity p xx is calculated as p xx = E x /J. 

linear MR. We have established a clear link between the stochas- 
tic dynamics of the electron trajectory around low- mobility islands 
caused by multiple scattering at the islands and the classical lin- 
ear MR in devices with different geometry and size. Although our 
observations have focused on the high-mobility inhomogeneous 
n-type InAs and In(AsN) epilayers grown on GaAs, our main find- 
ings and conclusions are relevant to a wide variety of structures, 
device geometries and material systems. Our results are therefore 
of general interest and will help guiding future research into the 
manipulation of the linear MR by tuning of the electron trajectories 
in tailor-made inhomogeneous materials and device structures. 

Methods 

Growth and fabrication. Our InAs epilayer (thickness t= 1. 5 jam) was grown by 
MBE on a semi- insulating (100) -oriented GaAs substrate. An In(AsN) epilayer 
(t= 1.3 um) containing a small concentration (N= 1%) of substitutional N-atoms was 
also grown by MBE using a radio frequency (RF) plasma N- source. The surface 
reconstruction was monitored by in situ reflection high-energy electron diffraction, 
while the substrate temperature was measured using an infrared pyrometer 
calibrated with the surface reconstruction transitions at a fixed As -flux. A Veeco 
UNI-Bulb RF plasma source was used as nitrogen source. It was set at a RF power 
of 160 W with a nitrogen flux of 5x 10 ~ 6 mbar for all epilayers. The N-content in 
the In(AsN) epilayers and the crystal quality of the layers were optimized by setting 
the growth temperature, T G , and the growth rate, r G , to given values in the ranges 
Tq = 400-460 °C and Tq = 0.3- 1 \im h~ 1 , and by keeping a minimum As-flux. The 
InAs epilayer was grown at T G = 460 °C. High-resolution X-ray diffraction measure- 
ments using a Bede QC200 diffractometer confirmed a relaxation of the crystal 
lattice by 98% in the InAs and In(AsN) epilayers. 

The wafers were spun with photoresist and patterned by photolithography 
into Hall-bars of length L= 1,250 um and width W= 250 jam. Two-terminal planar 
devices with InAs channels of length L = 2, 5 and 10 um and of width W= 1, 2 and 
5 \xm were also fabricated by electron beam lithography. Each channel was formed 
by dry-etching the epilayer down to the GaAs substrate. Metal contacts consist- 
ing of lOnm of Ti followed by 200 nm of Au were deposited onto the samples to 
provide Ti-Au ohmic contacts. 



Transport studies. Measurements of the current-voltage characteristics were 
performed in both direct current and pulsed mode. For the pulsed mode of 
acquisition of the current, a voltage pulse generator was used to generate saw- 
tooth single voltage pulses of duration of lOjis, and a digital oscilloscope was 
used for the synchronized and fast temporal recording of the applied bias and 
current, thus reducing excessive sample heating at high-operating currents. The 
current-voltage curves in pulsed magnetic field were taken within 1 ms at the top 
of the magnetic field pulse, during which the magnetic field remains constant to 
within ±0.5%. 

References 

1. Xu, R. et al. Large magnetoresistance in non-magnetic silver chalcogenides. 
Nature 390, 57-60(1997). 

2. Hu, J. & Rosenbaum, T. F. Classical and quantum routes to linear 
magnetoresistance. Nat. Mater. 7, 697-700 (2008). 

3. Branford, W. R. et al. Geometric manipulation of the high-field linear 
magnetoresistance in InSb epilayers on GaAs (001). Appl. Phys. Lett. 86, 
202116(2005). 

4. Friedman, A. L. et al. Quantum linear magnetoresistance in multilayer epitaxial 
graphene. Nano Lett. 10, 3962-3965 (2010). 

5. Zhang, W. et al. Topological aspect and quantum magnetoresistance of 
|3-Ag 2 Te. Phys. Rev. Lett. 106, 156808 (2011). 

6. Abrikosov, A. A. Quantum linear magnetoresistance; solution of an old 
mystery. /. Phys. A Math. Gen. 36, 9119-9131 (2003). 

7. Parish, M. M. & Littlewood, R B. Non-saturating magnetoresistance in heavily 
disordered semiconductors. Nature 426, 162-165 (2003). 

8. Parish, M. M. & Littlewood, P. B. Classical magnetotransport of 
inhomogeneous conductors. Phys. Rev. B 72, 094417 (2005). 

9. Demo, M. P., Yamamoto, S., Kasai, S., Ono, T. & Kobayashi, K. Large positive 
magnetoresistive effect in silicon induced by the space-charge effect. Nature 
457, 1112-1115 (2009). 

10. Sun, Z. G., Mizuguchi, M., Manago, T. & Akinaga, H. Magnetic-field- 
controllable avalanche breakdown and giant magnetoresistive effects in gold/ 
semi-insulating- GaAs Schottky diode. Appl. Phys. Lett. 85, 5643-5645 (2004). 

11. Heremans, J. Solid state magnetic field sensors and applications. /. Appl. Phys. 
26, 1149-1168 (1993). 

12. Patane, A. et al. Effect of low nitrogen concentrations on the electronic 
properties of InAsx_ x N x . Phys. Rev. B 80, 115207 (2009). 

13. Fahy, S., Lindsay, A., Ouerdane, H. & O'Reilly, E. P. Alloy scattering of n-type 
carriers in GaN^As^. Phys. Rev. B 74, 035203 (2006). 

14. Drachenko, O. et al. Cyclotron resonance mass and Fermi energy pinning in 
the In(AsN) alloy. Appl. Phys. Lett. 98, 162109-162111 (2011). 

15. Makarovsky, O. et al. Hot electron transport and impact ionization in the 
narrow energy gap InAs x _ X N X alloy. Appl. Phys. Lett. 96, 0521 15-0521 18 
(2010). 

16. Jacoboni, C. & Reggiani, L. The Monte Carlo method for the solution of charge 
transport in semiconductors with applications to covalent materials. Rev. Mod. 
Phys. 55, 645-705 (1983). 

17. Hamaguchi, C. Basic Semiconductor Physics 46 ( Springer- Verlag, London, 
2006). 

18. Fischetti, M. V. Monte Carlo simulation of transport in technologically 
significant semiconductors of the diamond and zinc-blende structures. 

I. Homogeneous transport. IEEE Trans. Electron Dev. 38, 634-649 (1991). 

Acknowledgements 

This work was supported by the Royal Society (UK), DFG KO 3743/1-1, AOBJ: 550341. 
We acknowledge Dr R. Airey (EPSRC, III-V National Facility, Sheffield) for processing 
our samples. 

Author contributions 

A.P., N.V.K. and O.M. designed and performed the experiments; N.M. performed the 
Monte Carlo simulations; L.E. took part in the discussions of the Monte Carlo studies; 
A.K. and Q.D.Z. grew the samples; A.P. wrote the paper; and all authors took part in 
the discussion/analysis of the data. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

Reprints and permission information is available online at http://npg.nature.com/ 
reprintsandpermissions/ 

How to cite this article: Kozlova, N.V. et al Linear magnetoresistance due to multiple- 
electron scattering by low- mobility islands in an inhomogeneous conductor. 
Nat. Commun. 3:1097 doi: 10.1038/ncomms2106 (2012). 

License: This work is licensed under a Creative Commons Attribution-NonCommercial- 
Share Alike 3.0 Unported License. To view a copy of this license, visit 
http://creativecommons.Org/licenses/by-nc-sa/3.0/ 



NATURE COMMUNICATIONS | 3:1097 | DOI: 10.1038/ncomms2106 | www.nature.com/naturecommunications 5 

© 2012 Macmillan Publishers Limited. All rights reserved. 



